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Abstract. We study a gas of TV hard disks in a box with semi- 
periodic boundary conditions. The unperturbed gas is hyperbolic 
and ergodic (these facts are proved for N — 2 and expected to be 
true for all N > 2). We study various perturbations by twisting 
the outgoing velocity at collisions with the walls. We show that 
the dynamics tends to collapse to various stable regimes, however 
we define the perturbations and however small they are. 
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Gas of hard balls is a classical model of statistical mechan- 
ics. The chaotic motion of its molecules can explain many 
phenomena observed in real gases and fluids. We study a 
2D gas of hard disks placed in a channel whose walls act as 
"Maxwell demon" that makes the angle of reflection different 
from the angle of incidence (twisting the outgoing velocity 
vectors). We present quite unexpected patterns developed by 
the gas dynamics under various twisting rules. 

1. Introduction 

The gas of hard balls is a classical model of statistical mechanics. 
Hard ball collisions produce strong scattering effect which makes the 
system behave chaotically and quickly relax to equilibrium (at least lo- 
cally). In mathematical terms, the gas of hard balls in finite container 
is widely regarded as a hyperbolic dynamical system (i.e., its Lyapunov 
exponents should not be zero) and its natural invariant measure (Li- 
ouville measure) is expected to be ergodic and mixing. 

More precisely, the celebrated Boltzmann-Sinai ergodic hypothesis 
states that the gas of TV hard balls in a container with periodic bound- 
ary conditions (a torus) is hyperbolic, ergodic, mixing (and Bernoulli). 
Attempts to prove this conjecture have long history [2EJ EHJ ED], and 
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at present it is proven under various conditions [13 [2U [23 [23 EZ], but 
not yet in its full generality. 

We note that the gas of TV hard balls on a d-dimensional torus has 
d+1 integrals of motion: its kinetic energy and its total momentum (a 
d- vector) are preserved, hence the position of the center of mass can be 
fixed. This leads to a great reduction of the phase space eliminating a 
total of 2d + 1 dimensions. 

If the container has rigid walls, the total momentum is no longer 
preserved, thus the phase space has higher dimensionality and is more 
complicated; such gases are even harder to study. Hyperbolicity and 
ergodicity have been proven only for N = 2 balls in a rectangular box 
(when the disks are no too large) [23] and for N > 2 disks in a very 
special 2D container with curved walls where each disk is confined to 
its own cell [5]. 

We are interested in a gas of N hard disks in a 2D rectangular con- 
tainer with partial periodicity, i.e., where two opposite walls are rigid 
but in the other direction boundary conditions are periodic (such a 
container can be regarded as a cylinder, rather than a torus). The 
hyperbolicity and ergodicity for such a gas are proved only for N = 2, 
see [22], but these properties are undoubtedly valid for all N > 2. 

The dynamics of this gas becomes more intriguing if small driving 
forces are added at the rigid walls, i.e., when collisions of the disks with 
the walls are modified by stochastic or deterministic perturbations. We 
consider deterministic perturbations where the angle of reflection is 
no longer equal to the angle of incidence, i.e., the velocity vector of 
the colliding disk is "twisted" a little right after the collision. Such 
twisting collision rules may appear physically unrealistic (they belong 
to the "Maxwell-demon" type of external forces), but they produce 
very realistic and interesting nonequilibrium phenomena such as shear 
flow and entropy production [7J. 

It is commonly expected (and observed empirically [7]) that the gas 
of hard balls under small perturbations remains chaotic and has a 
unique nonequilibrium stationary state, perhaps in the form of a Sinai- 
Ruelle-Bowen (SRB) measure, i.e., an ergodic measure with smooth 
conditional densities on unstable manifolds. These facts have been ac- 
tually proven for N = 1 particle in the the 2D periodic Lorentz gas 
under two types of small perturbations: external fields [3 QUI H2] and 
twisting collision rules [32]. Similar results were obtained for various 
classes of bounded billiard tables with twisting walls [U [T7] . 

We note that the unperturbed Lorentz gas with N = 1 particle is 
a uniformly hyperbolic dynamical system. Under small perturbations, 
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uniform hyperbolicity usually survives and makes the construction of 
an SRB measure possible. 

The gas of iV > 2 hard balls, on the other hand, is never uniformly 
hyperbolic. Non-uniform hyperbolicity can be easily destroyed even by 
arbitrary small perturbations making it hard to control the perturbed 
dynamics. For this reason there are no theoretical proofs of hyper- 
bolicity or the existence of SRB measures under any perturbations, to 
our knowledge. Moreover, one may expect that the perturbed gas is 
not always fully hyperbolic or ergodic, i.e., that there may be elliptic 
islands or multiple ergodic components, etc. 

The modest purpose of this paper is to point out that things are 
actually much worse, even for the gas of N = 2 disks. We show that if 
the disks are not too big (if their diameters are less than 0.5), then an 
arbitrarily small twist added to collisions at the walls tends to destroy 
the chaotic behavior of the gas causing a complete collapse of the dy- 
namics in a way that almost every phase trajectory converges to some 
trivial stable regimes. The collapse seems to happen for every type of 
small twists, though the limit stable regimes may be very different and 
sometimes quite bizarre. 

Most of the stable regimes we observed for N = 2 also appear for 
any number N > 2 of disks (provided their diameter is small enough). 
On the other hand, when the disks are not too small, then our stable 
regimes seem to disappear. We believe that if the diameters of the disks 
are greater than 1/N, then the gas is fully chaotic and SRB measures 
exist, in agreement with numerical evidence [7]. 

This work is motivated by discussions with J. Lebowitz and Ya. 
Sinai. Our original purpose was to prove hyperbolicity and construct 
SRB measures for at least some types of twisting collisions. However, 
we discovered that at sufficiently low densities all twists cause a total 
collapse (often to our surprise). Thus in the end we decided to report 
our findings in this paper. 

2. Model 

Our model is a system of N hard disks in a unit square 

D = {0 < x < 1, < ?/ < 1} 

that has rigid (reflecting) walls at y = and y — 1 and periodic 
boundary conditions at x = and x = 1. The disks are identical (have 
the same mass and radius), and they collide elastically with each other. 

We can represent our model as an infinite chain of copies of the 
square D placed in the infinite strip (channel) I = {0 < y < 1} where 
hard disks appear periodically; see Fig. [TJ 
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Figure 1. A periodic system of disks moving in a channel. 



We denote by qj = (xj, yi) the position of the ith disk and by pi = 
(ui, Vi) its velocity vector, i = 1, . . . , N. Now suppose a disk collides 
with a wall (either y = or y = 1). By the classical law its velocity 
P = Pbefore = (u,v) changes to Paftcr = (u, — v), and its kinetic energy 
i||p|| 2 is preserved, i.e., 

(2-1) 1 1 Pbefore 1 1 = 1 1 Pafter 1 1 - 

The total momentum p to tai = Pi + • • • + Pn of the system is not pre- 
served, but its x-component is preserved, i.e., the sum ui + ■ ■ ■ + 
remains constant in time. 

The entire (macrocanonical) phase space is 4iV-dimensional, and the 
energy surface is (47V— l)-dimensional. The conservation of U\ + - ■ -+un 
allows us to set it to zero, after which the x-coordinate of the center of 
mass will be constant and can be fixed, too. Thus the reduced phase 
space is (4iV — 3)-dimensional. The resulting system is expected to be 
hyperbolic (with 2iV — 2 positive Lyapunov exponents and the same 
number of negative ones), ergodic, and mixing. These facts have been 
proved for N = 2 in [23] . 

We now modify the law of collisions with the walls. We will preserve 
the kinetic energy of the disks, so that (12.11) still holds true. Thus it is 
enough to specify the angle of reflection if;, as a function of the angle 
of incidence, <p: 

(2.2) i> = f{<p). 

We measure the angles ip and if as shown in Figure [21 so that the range 
of our angles is the interval [0, n}. The function / does not depend on 
the point of collision, but it may be different for the bottom wall at 
y = and the top wall at y = 1. We will denote those two functions 
by f and fi, respectively. 
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Figure 2. Angle of incidence <p and angle of reflection ip. 

For classical collisions ip = <p, hence / is the identity function, /(</?) = 
tp. We will consider its small perturbations, i.e., our functions / will 
satisfy f(<p) = <p + small. For convenience we also assume that 

(a) /(</?) is continuous and strictly monotonically increasing; 



i.e., / is an orientation preserving homeomorphism of the interval [0, n]. 
This ensures that our dynamics will be invertible, i.e., every phase point 
has a unique past. 

As before, the energy surface in the phase space is (AN— l)-dimensional. 
Since our collision rules are translation invariant (independent of the 
x coordinate of the collision point), the energy surface is naturally fo- 
liated by invariant hypersurfaces on which the dynamics are identical. 
Thus we can factor this foliation out by replacing the x-coordinates 
of the disks by their relative x-coordinates x\ — x n , . . ., x n _i — x n , 
which eliminates one more variable and makes the (factored) phase 
space (4N — 2)-dimensional. Most of the time we will deal with N = 2 
disks, in which case the phase space will be 6-dimensional. 

Many collision rules (I2.2p quickly lead to trivial degenerate regimes. 
For example, if fk(p) < P for k = 0, 1 and all < <p < n, then 
^after = cos ip > cos <p = ^before, hence the re-component of the total 
momentum, u tota i = u\ + ■ — h will grow at every collision with the 
walls. It is then clear that all the particles will eventually move almost 
horizontally to the right (will be "blown away by wind"). 

To prevent the wind from blowing, we will restrict our study to 
collision rules where one wall counterbalances the effect of the other: 



(b) /(0) = and /(tt) = vr, 



(2.3) 



- /o(7T - p), 



i.e., the function fi at the top wall y = 1 acts exactly opposite to the 
function f at the bottom wall y = 0. Under this condition there is a 
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natural symmetry in the channel, hence drift in either direction (left 
or right) cannot be dominant in the whole phase space. 

3. One particle case 

To clarify the effect of our twisting collisions, we begin with the 
simplest case of one hard disk, i.e., we set AT = 1. Then the radius of 
the disk is irrelevant, and we can just make it a point particle bouncing 
between the two walls. 

Let (fo be the angle of incidence at the initial collision at, say, the 
bottom wall y = 0. Then the reflection angle ip = /o(<£>o) becomes the 
incidence angle <pi at the next collision at the top wall, i.e., <p\ = fo(<fo). 
Similarly, the incidence angle at the following collision at the bottom 
wall will be 

<P2 = fl((pi) = /l(/o(Vo)), 

and then the process will repeat periodically. By induction, 
V?2n = g n (<Po) for all n > 1, 

where g = f\ o / , thus the evolution of incidence angles is completely 
described by the iterations of the function g. For functions with oppo- 
site orientations, i.e., those obeying (12. 3p . we have 

(3-1) <7(v)=^-/o(*--/o(¥0). 

Note that g, just like fo, is an orientation-preserving homeomorphism 
of the interval [0,7r]. But not any homeomorphism of [0, ir] can be 
defined by (13. ip . so we describe the class of functions g satisfying (13. ip . 
To this end we introduce an involution j : [0, n] — > [0, it] defined by 
= ir — ip and note that 

9 = U° fo) (j ° fo) 

Thus, g = h 2 , where h = j o f is an orientation reversing homeomor- 
phism of the interval [0, 7r]. 

The map h: [0, ir] — > [0,7r] obviously has a unique fixed point <po £ 
(0, 7r), which automatically is a fixed point for g. For any other point 
<p ^ y?o there are exactly three possibilities: 

(a) h 2 ((p) = ip, then ip is a 2-periodic point for h and a fixed point 
for g; 

(b) \h 2 (ip) —ipo\ < \p — po\, then the images of <p under g will move 
toward tpo; 

(c) \h 2 (<p) —po\ > \<p — <po\, then the images of <p under g will move 
away from p ; 
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The non-fixed points of g make an open set, which is a union of disjoint 
intervals, we denote them by {I m }- The endpoints of each interval I m 
are fixed points for g, i.e., 2-periodic points for h (unless one of them 
is ip , of course). In each interval I m C [0, ir] all the points move under 
g in one direction - either toward (fio or away from (p . 

For each interval I m its image I m > = h(I m ) is another interval whose 
points move in the same direction (either toward p>o or away from ipo) 
as the points of I m . Note that J m = h(I m i) as well. We call the intervals 
I m and I m i dual to each other. Thus the intervals {I m } come in dual 
pairs. Dual intervals lie on the opposite sides of ipo, and all the points 
in both intervals move either toward p>o or away from ip . If I m , I m > 
and I n , I n i are two pairs of dual intervals, then one pair lies inside the 
other; see Figure [3) 



vim s 
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Figure 3. Dual intervals about the fixed point (fQ. 

If the number of intervals I m is finite, then it is necessarily even, and 
exactly half of them lies to the left of cpo, and the other half - to the 
right. The above structure is essentially a complete description of all 
interval maps g satisfying (13. lft : 

Proposition 3.1. Let g be an orientation preserving homeomorphism 
of the interval [0, ir] with an odd number of fixed points 

= y?_ fc < (f-k+i <••■<<#)<■-• <Pk-l < ¥k = 7T. 

Suppose that for each < % < k either all the points of both intervals 
(</?_j_i, tf-i) and (</?j, tfi+i) move under g toward (fio or all the points in 
these two intervals move away from ip . Then there exists an orienta- 
tion reversing homeomorphism h: [0, 7r] — > [0,vr] such that g = h 2 . The 
middle fixed point tpo is the (only) fixed point of h. 

The proof uses standard methods of one-dimensional topological dy- 
namics, and we omit it. 
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The central fixed point ip plays a special role, we will call it the 
center (of the map g); it is the only fixed point of h. Note that 

(3.2) fo(<Po) = - <^o, 

hence if a particle hits a wall at the angle ipo, it turns around and flies 
straight back. Its trajectory is then periodic not only in the angular 
coordinates, but also in the spatial coordinates; see Fig. HJ 




Figure 4. A periodic trajectory running between two collisions. 

4. A SPECIAL FAMILY OF COLLISIONS WITH A TWIST 

While most of our conclusions apply to generic functions /o and fi 
satisfying (12.3| . it will be convenient to use one particular family of 
functions to clarify our arguments. 

First, it is convenient when the center (po is at the geometric center of 
the interval [0, tt], i.e., <po = tt/2. Then the special periodic trajectories 
described in the end of Section [3] move vertically, up and down, i.e., 
they just bounce between the walls like the regular periodic billiard 
trajectories. 

Next, for simplicity we consider functions g that have one fixed point 
<Po — 7r /2 (other than and tt), which will automatically be the center. 
Then there are three fixed points total: 0, tt/2, and tt. The intervals 
(0, 7r/2) and (tt/2, tt) either both move toward tt/2 or both move away 
from 7r/2. In the former case tt/2 will be a stable fixed point and attract 
the entire interval (0,7r). In the latter case tt/2 will be an unstable 
(repulsive) fixed point and then each interval (0, tt/2) and (tt/2, tt) will 
be attracted by its other endpoint, i.e., by or tt, respectively. We 
will say that in the former case the map g has a stable center and 
in the latter - an unstable center. These cases are quite different, in 
dynamical terms, and both are interesting. 

We note that in order to make tt/2 a fixed point for g satisfying (13. ip . 
we need to set f^ir/2) = tt/2, according to (13.21) . Then tt/2 is a fixed 
point for fo as well. 

Thus fo must have three fixed points: 0, tt/2, and tt. It is also 
convenient if the time reversal of the collision rule defined by /o belongs 
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to the same family of collision rules. The time reversal collision rule is 
defined by <p \-> /^(v 9 ) & t the bottom wall and (p \-> fi(<f) at the top 
wall, where the functions and /f satisfy 

(4-1) fo=f{ 1 and fT = fo l 

in accordance with (12. 3p . 

In view of the above requirement we can define fo by 

(4.2) tan/ (</>) = e A tany?. 

Then, due to (jOjl . 

(4-3) fM = Mv>) 

so both walls obey the same collision rule! Note that there are no 
restriction on A here, our rules work well for any A G (— oo, oo), but for 
the twist to be small we will only consider A ~ 0. 
Due to fl4.ll) . the time reversal collision rules satisfy 

tan/o~(<£>) = tan f\{f) = e _A tany9, 

hence it belongs to the same family of rules (I4.2|) . but A must be re- 
placed with —A, i.e., reversing the time corresponds to negating A. 
A direct differentiation of (14. 2[) gives 

ti( m \ - „-a sin 2 / (y) _ \ cos 2 / (y) 
sin Lp cos^ if 

Hence 

= e~ x and /^(0) = ^(tt) = e\ 

Thus, A > corresponds to a stable center and A < to an unstable 
center, i.e., we have dynamics of both types. 

We note that the rules (I4.2p - (j4.3p can be rewritten as follows: at 
every collision with the wall the incoming velocity (u~,v~) and the 
outgoing velocity (u + ,v + ) are related by 

(4.4) u + /\v + \=e- x u-/\v-\. 

This makes it convenient for numerical simulations: one can recompute 
the velocity vectors by (14. 4p without using the angle cp or its tangent. 

5. Collisions with a stable center 

We begin with a single disk and a small positive A > 0. Again, we 
assume that the disk has zero radius, i.e., it is just a point particle. 

Recall that /q(7t/2) = e _A < 1. Thus if the angle of incidence is close 
to 7r/2, i.e., if = 7r/2 + 5 for some small 5, then after n collisions at the 
walls it will be <p n = ir/2 + S n , where S n — > 0. 
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More precisely, as the disk moves between collisions from wall to wall, 
its displacement in the horizontal direction is given by Ax = cot tp = 
tan 6. Due to fl4.2p - fl4.3p . after the next collision its displacement will 
be e~ x cot </? = e~ x Ax. Thus the displacement in the x direction is 
precisely decreasing by a factor e~ x < 1. 

Thus the displacements in the x direction make a geometric progres- 
sion and the total displacement is 



The particle's trajectory converges to a vertical line and the velocity 
vector aligns vertically at an exponential rate. 

Consider now the system of N disks. We assume that the diameter 
d of the disks is small enough so that the disks can be lined up along 
one vertical or horizontal line in the unit square, i.e., d < 1/N. 

For any disk with center at (xi, yi) and velocity vector (m, Vi) denote 
by Si the total displacement in the x direction of that disk if it is allowed 
to bounce between the walls alone, as if the other disks did not exist. 
By (15. ip . Si is finite, and in fact 



Let U be the the subset in phase space satisfying the following con- 
dition: 

(5.3) \xi — Xj | > d + Si + Sj Wi ^ j 

(note that the x coordinates must be taken modulo 1, as we have 
periodic boundary conditions at x = and x — 1). 

Now (15 .3p guarantees that the projections of the disks onto the x-axis 
do not overlap and will not overlap at any time in the future. Thus our 
N disks will never collide with each other, all of them will be bouncing 
between the walls and the trajectory of each disk will converge to some 
vertical line. 

In this limit regime, the x-component of the velocity vector of each 
disk is zero. Thus those limit trajectories form a family of codimension 
N in phase space. Note that the phase space of our system is (AN — 1)- 
dimensional (the only constraint is the conservation of the total kinetic 
energy). Thus our family of stable limit trajectories is a (3N — 1)- 
dimensional submanifold in the phase space. We denote it by S. 

Every phase trajectory starting in U converges to S at an exponential 
rate. The set U is invariant, i.e., all the trajectories originating in U 
stay there in the future. Note that U is an open subset of the phase 



(5.1) 
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Figure 5. Convergence the a stable limit regime. 

space, thus it has a positive Lebesgue measure. Therefore the limit set 
S has an open basin of attraction of positive Lebesgue measure. 

It is more convenient to work with the collision space Ai which con- 
sists of all phase states such that either some two disks collide or a disk 
collides with a wall. (At the moment of collision, the velocity vectors 
change discontinuously, and it is customary to include in Ai only the 
postcollisional velocity vectors). The induced map T — T\: Ai — > Ai 
is called the collision map. 

We denote by vq the normalized measure on Ai that is invariant 
under the map T corresponding to A = 0. This map corresponds to 
the classical specular reflections at the walls (where ip = <p), so T 
is the collision map for the classical gas of hard disks at equilibrium. 
The measure vq is absolutely continuous with respect to the Lebesgue 
measure on Ai, and it has a strictly positive density. At the collisions 
with the walls, the density is proportional to sin <p [8] . 

Let Uq = U fl Ai denote the subset of the collision space where 
( 15. 3p holds. It is invariant under T, i.e., T(U ) C U . Its trajectories 
converge to the set So = S fl Ai, i.e., U is the basin of attraction of 
Sq under the map T. The basin of attraction Uq may be regarded as 
a trap, or a "hole" - every phase trajectory that enters it will never 
come back. 

6. Holes and escape rates 

Dynamical systems with holes have been studied extensively in the 
past decades, both numerically and theoretically [2j |9j HU |20j [22] • In 
particular, many researchers studied billiards with holes [31 HI [131 ESJ 
|19] . It is generally observed (and in many cases proved) that holes at- 
tract almost the entire phase space, i.e., almost every trajectory sooner 
or later enters the hole and never returns (in other words, it escapes). 
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This phenomenon is also interpreted as the leakage of mass (phase 
volume) through holes so that the remaining phase space gets thinner. 

The holes are usually characterized by the escape rate p > 0, which 
basically says that the fraction of phase space that has not escaped 
through the holes (not leaked out) by the time n > (where n is 
discrete time, i.e., the collision counter) decays exponentially at a rate 
p, i.e., it is of order e~ pn . More precisely, if M. n denotes the subset of 
the phase space Ai that has not escaped through the holes by the time 
n, then 

lim - lnLeb(A^ n ) = — p. 

n— >oo n 

For any phase state X G M. we denote by t(X) the escape time, 
i.e., the time the trajectory of X enters the hole. Then t(X) has an 
approximate exponential distribution on M. with parameter p (with 
respect to the Lebesgue measure), and its average is 1/p. 

Many studies investigate how the escape through holes is affected 
by the size of the latter. When holes get smaller (shrink), the escape 
rate decreases and the escape time grows. In our case, the "hole" Uq 
described above depends on the parameter A > 0. When A decreases, 
the hole Uq shrinks and in the limit A — > the hole converges to the 
sub manifold So C M. of zero volume. 

We have investigated the escape process in our model by numerical 
simulations. We used N = 2 disks of diameter d = 0.1 with the collision 
rules at the walls defined by (I4.2p - fl4.3l) for various small A > 0. We 
have estimated the "escape time" numerically, for different A's. For 
each A we simulated the dynamics from 10 5 randomly chosen initial 
states, stopping the simulations whenever the system escaped into the 
"hole" U , and computed the mean escape time p T . 

It naturally grows, as A gets smaller, and Figure shows how its 
grows on the log-log scale. The plot clearly demonstrates a linear 
pattern suggesting that p T ~ \~ b for some b > 0. The least squares 
fitting line in Figure[6]has slope 2.11, so we may guess that b « 2. Below 
we give a heuristic argument supporting (and refining) this conjecture. 

Recall that the "hole" Uq C M. is an open set of positive Lebesgue 
measure, and it is invariant under the collision map T, i.e., T(U ) C Uq. 
Hence the sets 

U x = T-\Uq) \ Uq and U k =T~ k+l {U 1 ) 

for k = 2,3,... are disjoint and each of them has positive Lebesgue 
measure, too. Their union Uoo = U^ £4 contains all the phase points 
I 6 M that end up in the hole Uq eventually. Our first conjecture 
(well supported by our numerical observations) is that almost every 
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Figure 6. Plot of In [i T versus — In A. The least squares 
line has an estimated slope of 2.11. 

phase point X G M. eventually ends up in the hole, i.e., Uoo covers the 
entire collision space M. (up to a null set), i.e., M. = (mod 0). By 
direct analysis (we omit details) 

(6.1) Leb(f/ ) = 0{X N ) and Leb(tfi) = 0{X N+1 ). 

and similar estimates hold for vq(Uo) and vq(U\). 

Next, the measure z/ is invariant under the original, twist-free col- 
lision map Tq, but not for A ^ 0. In fact, it is preserved by the inter- 
particle collisions, but at collisions with the walls it gets compressed 
if f' (tp) < 1 and gets expanded if /o(y?) > 1. In our case, /q(v) < 1 
whenever c\ < <p < c 2 for some constants < c\ < c 2 < n whose values 
are not essential. 

For any point XeMwe denote the Jacobian of the inverse collision 
map T" 1 , with respect to the canonical measure v$ by J(X) = e" / ^ x \ 
If X corresponds to an inter-particle collision, then J(X) = 1, hence 
j(X) = 0. Otherwise 7 has a small (of order A) value. 

Now for every point X G U\ and k > 1 we denote 

J k (X) = J{T~ k+ \X)) = e^, lk = 7 (T^ +1 (X)), 
then we have 

(6.2) MU k+ i)= [ J 1 ---J h du Q = f e* + -+* du Q . 
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Next note that the sequence j± , . . . , 7^ is determined by the colli- 
sions of the disks with the walls along the past trajectories of points 
X G U\. As time runs backwards, the disks begin colliding with each 
other and a chaotic regime quickly sets in. Then the distribution of 
every trajectory in the phase space will be fairly close to uniform (equi- 
librium). Thus the values of 7 n could be treated as nearly independent 
random variables. They are all of order A, hence the sum of k of them 
can be expected to grow as Xy/k, in the spirit of the the Central Limit 
Theorem. Thus by (I6.2p we can expect that 



v (U k+1 )~ / e AVfe dX~A JV+ V 

JUi 

as long as the chaotic regime continues. The obvious limitation fo(f/fc+i) < 
1 gives us an upper bound for k: 

, x , (In A) 2 

(6.3) k < const v ■ / 

X 2 

Perhaps it is more reasonable to expect that z/ (?7 fc+1 ) ~ 1/k, but this 
would give us pretty much the same upper bound (16. 3p on k. 

Thus we see that the chaotic regime, when the particles collide 
and the measure of the regions U k tends to grow, lasts for about 
£}((ln A) 2 /A 2 ) collisions. As a result, the mean 'escape time' is of the 
same order: 

(6.4) /i T = (9((lnA) 2 /A 2 ). 

On the log- log scale x = — In A and y = In fi T adopted in Figure [6] this 
means 

(6.5) y = 2x + 2 In x + const. 

Accordingly, we used a functional relation y = ax + b In x + c to ap- 
proximate our data plotted in Figure El and the least squares fit gives 

y = 1.985x4- 0.233 lnx + 0.811. 

The estimated slope of 1.985 is in a good agreement with the theoreti- 
cally predicted slope of 2. The logarithmic coefficient 0.233 is not close 
to 2 in (16. 5p . but it is of secondary importance for the estimate (16. 4p . 

We now investigate the dynamics beyond the bound (16. 3 p . The mea- 
sure Ho(Uk) cannot increase with k forever; in fact we must have 

A*o(£4) — > as k 00, 

and moreover, fj,o(Uk) should decay fast because the series ^2 k f^o(Uk) 
is summable. 
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Thus for large k (those far exceeding the upper bound ( 16. 3ft ) the 
dynamics cannot be chaotic: it is dominated by the collisions with the 
walls that compress the phase volume (when the time runs backwards) . 
Those collisions have incidence angles close to or ir, so that the par- 
ticles hit the walls nearly tangentially. This means that the particles 
move nearly horizontally in the channel. 

We recall that running the time backward for a dynamics with pa- 
rameter A is equivalent to reversing the velocity vectors of the disks 
and running the time forward for the dynamics with parameter —A. 

This suggests that for the dynamics with negative parameter A < 
the limiting regime also exists, and it involves the disks flying al- 
most horizontally and experiencing nearly tangential collisions with the 
walls. This will be established, with some rigor, in the next sections. 

7. Collisions with an unstable center 

Here we analyze the dynamics of N disks with a negative parameter, 
A < 0. In this case ip = and ip = n are stable fixed points for both 
functions fo and f\ and the derivative at these points is 

/o(0) = /oW=e A <l 

Thus if the angle of incidence is close to or ir, i.e., ip = 5oTip = 7i — 5 
for some small 5, then after a collision at the wall it will be even closer 
to or ir, i.e., it will be (pi — 5\ or tpi = ir — 5\, respectively, with 
5\ < 5e x with some constant A' < (A' ~ A for small S). 

However the angles of incidence cannot converge to or ir monoton- 
ically in a manner similar to their convergence to ir jl in the previous 
section. Indeed, after a nearly tangential collision with the wall the 
disk has to move across the channel to the opposite wall, and now on 
its way across it will be very likely colliding with other disks. 

Lemma 7.1. With probability one the disks will continue colliding at 
arbitrary distant future. 

Proof. Suppose that the disks never collide with each other after some 
time T > 0. Then they just collide with the walls, and hence their ve- 
locity vectors align almost horizontally and converge to some horizontal 
vectors. 

By simple geometry, the disks can avoid each other on their way 
from one wall to the other only if the horizontal components of their 
velocities are nearly equal. And as the vertical components get smaller, 
the horizontal components would have to get also closer to each other. 
So in the limit, as time grows to infinity, the velocity vectors of the 
disks would have to converge to a common limit! 
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However, the collisions with the walls do not alter the kinetic energy 
of the disks. Thus if the disks somehow managed to avoid each other, 
the speed of each disk would remain constant. And since their velocity 
vectors must have a common limit, their speeds must have been equal 
all the time! In other words, the 'equal speeds' situation should have 
occurred right after the last collision between the disks in the past. 
This event is exceptional and occurs with probability zero. □ 

Thus almost surely our disks will keep colliding with each other for- 
ever. For this reason the limiting stable regime (if one exists) could 
not be as simple as the one we have seen in the previous section. Still 
a limiting stable regime exists, as we will show next. Again we assume 
that the disks are not too large. 

Let us set the total kinetic energy to iV/2, so that the velocity vectors 
(ui,Vi) of the disks will satisfy 

(7.1) ul + vf + --- + u 2 N + v 2 N = N. 
Now we define a subset W + in phase space by 

(7.2) W+ = {ui + ■ ■ ■ + u N > y/N(N-l)}. 
An elementary calculation shows that for every X G W + 

(7.3) Ui > for all i = l,...,N. 

Now consider the trajectory of a phase point X e W+. At every 
collision with a wall, the incoming velocity vector (v,i, Vj) points to the 
right, due to ( I7.3jl . Then the postcollisional velocity vector (u'^v^) will 
turn toward the wall, i.e., u\ > and \vl\ < \vi\. At every collision 
between disks i and j, the velocity vectors of both disks may change, 
but their total momentum is preserved, i.e., Ui + Uj remains the same. 
We see that the total horizontal momentum of the system 

(7.4) M u = Ul + ■ • ■ + u N 

is monotonically increasing: it is preserved at every interparticle colli- 
sion and increases at every collision with a wall. Thus the trajectory 
of every point X e W+ will remain in W + forever. Furthermore, the 
value of M u will grow and converge to a limit. 

We illustrate this process for N = 2. Figure [7] shows W + in the 
U\,U2 coordinates: W + is bounded by the circular arc u\ + u 2 = 2, 
Ui, U2 > 0, and its chord U\ +U2 = V%- Then W + is foliated by parallel 
lines u\ + U2 = c, a/2 < c < 2. The point (ui,u 2 ) stays on the same 
line U\ + «2 = c at every collision between the disks and moves up to 
another line u\ + u 2 = d with d > c during every collision with a wall. 
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Figure 7. The region W + and the dynamics inside it. 

Next we investigate the limit regime. Clearly there is a limit chord 
ui + u 2 = Coo < 2. If Coo = 2, the chord degenerates to a single point 
ui = u 2 = 1. In that case vi = v 2 = 0, so the limit regime consists of 
both disks moving horizontally with the same unit speed. 

If Cqc < 2, then the points (^1,^2) along the trajectory of X e W + 
accumulate on the chord u\ + u 2 = Cqo. We claim that they actually 
converge to one of the end points of that chord. Indeed, when the point 
(ui, u 2 ) is at a distance 5\ > from the end points of the current chord 
ui + u 2 = c, then u\ + u 2 < 2 — 5 2 for some 5 2 > (which depends on 
81). Therefore \vi\ > 5% or \v 2 \ > 5 3 for some 5 3 > (which depends 
on S 2 ). Now we claim that at the next collision with a wall by one of 
the disks its vertical velocity will be \vi\ > 64 for some S4 > (which 
depends on 5 3 ). Indeed, a sequence of consecutive collisions between 
the two disks can be reduced to a simple billiard trajectory in a 2D 
periodic Lorentz gas [Tlj Section 4.2], and then one can easily check 
that the vertical components of their velocities cannot both vanish at 
the time they hit the walls. 

Next, when a disk with a vertical velocity |t>j| > 64 hits a wall, its 
velocity vector will be rotated so that U{ will grow by some 65 > 
(which depends on 64). Thus the point (ui, u 2 ) will move up to a chord 
u-i + u 2 = c + 5q for some 5 6 > (which depends on 5 5 ). Hence the 
convergence to a limit chord can only occur when (ui, u 2 ) converges to 
one of its end points (and, as a result, v 1 , v 2 converge to zero). 

In the above limit regime, when Cqo < 2, both disks move horizontally 
but at different speeds, u\ 7^ u 2 . We believe this actually happens 
with probability zero, but we can only give a heuristic argument. By 
Lemma 17.11 the disks have to collide from time to time. When they 
collide, their relative horizontal velocity is u\ — u 2 ~ u\ — u 2 ^ 0, but 
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their relative vertical velocity is V\ — i>2 = s ~ 0. The collision of two 
disks can be reduced to a 2D periodic Lorentz gas with infinite horizon 
[IH Sect. 4.2]. 




Figure 8. Nearly horizontal trajectories in a periodic 
Lorentz gas. 

It is known in the studies of periodic Lorentz gases that if the vertical 
velocity of the moving particle before the collision is e, then after the 
collision it is typically of order y/e, i.e., v[ — v' 2 = 0(y/e); see Figure El 
More precisely, there are p, q > such that the relative measure of 
initial conditions for which v[ — v' 2 = 0(e 1 ^ 2 ~ p ) is 0(e q ); see [31]. Thus 
in the course of infinitely many successive collisions between the disks 
their vertical velocities would explode sooner or later with probability 
one, violating the assumption < 2. This argument is not quite 
formal, as the word "probability" here refers to the canonical measure 
b>o, and in our dynamics the images of this measure keep changing, but 
we believe the conclusion is correct. 

To summarize, we proved that in the limit regime the particles move 
horizontally to the right. We also conjecture that with probability one 
their horizontal velocities are equal: 

(7.5) u\ = ■ ■ ■ = un = 1 and v± = ■ ■ ■ = vn = 0. 

In other words, almost all trajectories in W + converge to a limit regime 
where all the particles move horizontally and at unit speed. This is a 
stable limit regime, it corresponds to a submanifold S + . This subman- 
ifold is defined by equations (I7.5p . hence its dimension is 2N (corre- 
sponding to the free coordinates (x^f/j) of all the particles). 
There is a symmetric region W_ in phase space defined by 

(7.6) W- = {ui + • • ■ + u N < —\JN(N - 1)}, 

where all the particles move in the negative x direction. By a similar 
argument, almost all trajectories X G VF_ converge to the stable limit 
regime where all the particles move horizontally to the left, and (we 
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again conjecture) at unit speed: 

(7.7) u\ — ■ ■ ■ — un — — 1 and V\ — ■ ■ ■ — Vn — 0. 

This corresponds to a submanifold SL of dimension 2N. 

Note that we have two stable limit regimes. Each one attracts a part 
of M. that has a positive Lebesgue measure. Their basins of attraction 
are obviously disjoint, and due to symmetry they must have the same 
^o-irieasure. Thus we have a coexistence of two attracting mechanisms 
in phase space. 



empirical data + 
0.983737 x+1 .21 81 74 
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Figure 9. Plot of ln/i T versus — ln|A|. The least 
squares fit has slope of 0.983737. 

Just like in Section we estimated the escape time numerically. We 
used N = 2 disks of diameter d — 0.1 with the collision rules at the 
walls defined by f l4.2p - fl4.3p for various small A < 0. Figure M shows 
the mean escape time /i T versus A on the log-log scale. The plot clearly 
demonstrates a linear pattern and the least square fitting line has slope 
0.984 suggesting that ji T ~ |A| _1 . Thus typical trajectories now escape 
much faster than in the case of stable center. We discuss the reason 
for the faster escape below. 

Our main argument in Section used the central limit theorem and 
was based on the chaotical character of the motion of the balls (before 
they enter the stable regime). However when the balls collide with the 
walls, our twisting rules fl4.2p - fl4.3p with A < make the horizontal 
components of their velocities increase, thus the x component of the 
total momentum, u\ + «2, tends to grow (in absolute value). This 
tendency causes a drift toward larger values of \u\ + M2I, i.e., toward 
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the trapping regions W±. Thus we should regard u\ + u 2 as a one- 
dimensional Ito diffusion process with a non-zero drift [T8"l I2T] . 



Average change of the x-velocity at collisions with the walls 
0.0025 I 1 1 1 1 1 1 1 1 




-0.0025 1 1 1 1 1 1 1 1 1 1 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

x-velocity 

Figure 10. Average change of the x component of the 
disks velocity u versus its current value. 

Figure (TU] described the degree of the drift: it shows the average 
change of u, the x component of the disk velocity, versus its current 
value (the plot was computed for A = —0.015). We see that the ab- 
solute value \u\ always tends to increase, and this tendency is strong 
everywhere except when u is close to or ±1. In Section [6] we dealt 
with the balls that were approaching the hole Uq; then they moved 
nearly vertically, hence the drift was small and could be ignored. Now 
the balls are approaching the regions W±; thus the drift is almost at 
its peak and its effect is crucial. Since the average drift, per collision, 
if of order |A|, typical trajectories reach the regions W± of order |A| _1 . 

Figure [ID] also shows that the drift persists as long as |u| < 1, though 
it decreases as \u\ — > 1. This supports our conjecture that typical 
trajectories converge to a limit regime where all the particles move at 
the same unit speed. 

8. Generalizations 

To summarize, the dynamics with unstable center (A < 0) is quite 
different from the one with stable center (A > 0). It has two distinct 
limit regimes (5+ and S-), both have dimension 2N, while in the other 
case we had a single stable regime S with dimension 3 N — 1 . 
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The attracting regions W± ("holes") are of fixed size and measure 
(independent of A), while for stable center the "hole" Uq was quite 
narrow and had measure 0(X N ). 

Next we generalize our results. Our analysis of the stable center case 
applies to any stable fixed point of the map g = f\of inside the interval 
(0, tt). Every such point produces a T-invariant open set Uq C M., i.e., 
T(Uq) C Uq, which acts as a "hole". It has a basin of attraction in Ai 
of positive Lebesgue measure, and the limiting stable regime consists 
of disks moving with periodic incidence angles (assuming that the disks 
are small enough to avoid colliding with each other). 

Of course if the map g : [0, tt] — > [0, tt] has more than one stable fixed 
point, the corresponding basins of attraction are disjoint, so we have a 
coexistence of several attracting regimes. 

Our analysis of the unstable center case applies to any map g with 
stable fixed end points, and tt. In that case we have two attracting 
regimes similar to S + and 5_ above, and each has its own basin of 
attraction. 

We must note that in the general case the "holes" W± need be defined 
more cautiously than (17.21) and (I7.6p . Precisely, we need to choose a 
small e > and define W± by 

W± = {±(«i H h Un) > N — e }. 

It is easy to verify that there is a c = c(eq) > such that for every 
phase state in W± we have 

±Ui > 1 — c and \vi\ < c 

for all % — 1, . . . , N. Moreover, c — > as e — > 0. So we can choose, 
for example, c = 0.1 and fix the corresponding e > in the definition 
of W± . Then of our analysis of the convergence to S± will not require 
significant changes. 

In summary, every stable fixed point of the map g : [0, tt] — > [0, tt] 
leads to a stable regime that attracts a set of positive Lebesgue measure 
in phase space. 

In our numerical experiments, we mostly observed stable regimes for 
N = 2 disks. Our analysis shows that they exist for larger iV's, too, but 
they become increasingly difficult to observe experimentally. Indeed, 
the sizes of the holes in phase space are exponentially small in iV (cf. 
( 16. ip . and a similar estimate can be derived for (17. 2p ). hence the escape 
time grows exponentially with N. Thus in physical systems with a large 
number of molecules the stable regimes become almost unobservable. 
Still, they exist for any N, and they dominate the dynamics for small 
iV's. 
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One may wonder whether attracting regimes exist when g has no 
stable fixed points. For example, what if g((p) = f is the identity map? 
The next section presents the most striking result: even in that case 
stable regimes may exist, which attract almost every phase trajectory! 

9. Time reversible collisions with walls 

Of special physical interest are twisting collisions (12. 2 p that make the 
dynamics time reversible. This means that if we reverse the velocity 
vectors of all our disks, then they will move along their past trajec- 
tories backwards. By direct inspection, the collision rule (I2.2p is time 
reversible if and only if the function / satisfies 

(9.1) f(ir-f(<p))=ir-<p. 

This condition implies that the graph of the function ip = f(ip) is 
symmetric about the line ip = tt — if. 

If our collision rule is time reversible, i.e., satisfies (19. ip . then we 
immediately arrive at g(tp) = f. 




Figure 11. Collision rule pushing the particles to the right. 

A simple example of a time-reversible rule is 

(9.2) cot f k (<p) = (-l) fc A + cot <p 

One can easily see that this definition is consistent with (12 .3p and (19. ip . 
The rule ( 19. 2ft can be written in the notation of ( 14.41) as 

(9.3) u + /\v + \ =u-/\v-\ + (-l) k X. 

This formula has a simple geometric interpretation: suppose a particle 
collides with a wall at a point with x-coordinate Xq, and we extend its 
trajectories before and after the collision until they cross the other wall 
at points whose x-coordinates we denote by X-\ and x\, respectively 
(see Figure [Til, then 

Xl - X = X - X-i + (-l) fc A. 
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If, for example, A > 0, then the above relation implies that the bottom 
wall y = pushes the particles to the right, and the top wall - to the 
left. This creates shear flow in the channel of the type studied in [TJ. 

Next we describe a special regime in the above shear flow for N = 2 
disks. Suppose the disks move with opposite velocity vectors (i.e., 
their total momentum is zero) and they collide with opposite walls 
simultaneously. Then after the collision they again move with opposite 
velocity vectors. 

When the disks collide with each other, then by symmetry their 
point of contact has coordinate y = 0.5, i.e., the collision occurs right 
in the center of the channel. After the collision the disks move with 
opposite velocity vectors again, and will collide with the opposite walls 
simultaneously, etc. Hence this regime is invariant under our dynamics. 
See illustration in Figure [T2l 



Figure 12. Special regime S* with symmetric motion. 

We denote the above special family of phase trajectories by S*. Note 
that dimS 1 * = 3, while the (reduced) phase space is 6-dimensional, cf. 
Section H 

A striking fact we discovered by numerical simulations is that the 
regime S* is stable and attracting. Almost every randomly selected 
phase trajectory eventually stabilizes near S* and then evolves in a 
symmetric fashion so that the disks move with opposite velocity vectors 
at equal distances from the opposite walls. 

Actually, the regime S* may be also stable for the dynamics with an 
unstable center (Section [7J, but this happens only for relatively large 
perturbations (A > 0.6). Since we are primarily interested in small 
perturbations, we will not discuss this last fact. 

We provide a semi-heuristic argument showing that the regime S* 
is stable for the time-reversible collision rules, i.e., phase trajectories 
near S* tend to get closer to S* in the future. 
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We perturb the symmetries of the regime S* and show that pertur- 
bations tend to decrease, on average. There are two symmetries in the 
regime S*: the velocity vectors of the disks are opposite, v and — v, 
and their y-coordinates sum up to one, i.e., y\ + y 2 = 1. 

First we perturb the velocity symmetry, i.e., suppose that the ve- 
locity vectors are v + £ v and — v, i.e., the total momentum is a small 
<5v ^ 0. The interparticle collisions do not change the total momentum. 
When the disks collide with (opposite) walls, then the postcollisional 
velocity vectors will be denoted by v' + 8v' and — v'. Our goal is to 
show that \5v'\ tends to be smaller than \8v\, on average. 

We decompose 5v = 8v» + 5vj_ into the components parallel and 
perpendicular to v, respectively. Similarly, 5V = 5vi + 5v' ± are the 
components of <5v' parallel and perpendicular to v'. Due to the conser- 
vation of the kinetic energy of each disk at the collision with the wall 
we have 

|v'| = |v| and |v' + 5v'\ = |v + 5v\ 

hence 

(9.4) |5v(|| = \5v\\\, 

to the leading order. We will show that \5v' ± \ < |5vj_|, on average. 

Let ip and (p+S(p denote the directional angles (i.e., angles made with 
the positive ) of the velocity vectors v and v + 5v, respectively. 

Let ip and ip + 5ip denote the directional angles of the velocity vectors 
v' and v' + 8v', respectively. Then 



to the leading order, hence 

(9-5) Ml = l^t = |/^)| 

\5v±\ \d(p\ 

(recall that fk is monotonically increasing, i.e., f' k > 0). 

Now after the collisions with the walls the particles move across the 
channel, and they either collide with each other or they miss each other 
and hit the opposite walls. In the latter case we denote by v" + 5v" and 
— v" their new postcollisional velocity vectors, and again decompose 
5V = 5v(f + 5\" ± . Then, inductively, 



and 



\5V[_\ 


\5v'i\ 




\8v±\ 


\5V ± \ 


\8v±\ 



|5v[[| = \8v'u | = \5v\\ | 
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therefore \5v"\ = \8v\. 

We see that as long as the particles collide with the walls only, the 
perturbation vector 5v changes periodically, with period two. So for 
an even number of wall collisions between two successive interparticle 
collisions, the net result will be zero change, i.e., \8v\ will remain the 
same. For an odd number of wall collisions between two successive 
interparticle collisions, the net result will be the same as for just one 
wall collision, i.e., (19. 5p will apply. Next we relate \8v±\ to |5v|. 

Consider an interparticle collision. It preserves the entire vector 8v, 
but the direction of the postcollisional velocity vector v can be regarded 
as a random variable uniformly distributed in the entire range [0,27r], 
due to the scattering nature of the elastic collisions of hard disks. Let 
j3 denote the directional angle of the vector 8v. Let tp denote, as 
usual, the directional angle of the outgoing velocity vector v. Then 
\8v±\ = \8v\ | sin(/3 — tp)\. The perturbation 8v' at the next interparticle 
collision will be 



where x = [f' k ] 2 for an odd number of intermediate wall collisions and 
x = 1 for an even number of intermediate wall collisions. 

Next we estimate the total change of the norm \8v\ over a long period 
of time (0, T) during which n interparticle collisions occur. We have 
n — 1 intervals between successive interparticle collisions, and some of 
them (say, m < n — 1 of them) have an odd number of collisions of each 
particle with the walls, while others (i.e., n — m — 1 intervals) have an 
even number of collisions of each particle with the walls. Our previous 
analysis can be summarized as 



(9.7) log|5v T |-log|<5v | = \ ^log[l-(l-[4(^)] 2 )sin 2 (A-^)], 



where the summation is taken over the intervals with odd numbers of 
wall collisions. 

Due to the randomization caused by the scattering effect of the in- 
terparticle collisions we can treat the angles 0i and </?j's as independent 
random variables with uniform distribution in their ranges < fli < 2n 
and < tfi < ir. We will prove in Appendix that the average value of 
each term in (19. 7p is negative: 

Lemma 9.1. For both k = 0, 1 we have 



(9.6) 



\8v'\ 2 = |5v| 2 (cos 2 (/3 -<p) + xsin 2 (/3 - <p)) 



in, 



i=l 




log[l - (1 - [f'M?) sm 2 (/3 - ip)] dpd V = ^ < 0. 
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Thus the sum in ( 19. 7p approaches — oo linearly in m (i.e., linearly in 
time), hence the norm of the perturbation <5v decreases exponentially 
in time. 

Second, we perturb the other symmetry of the regime S*, i.e., we 
suppose that an interparticle collision occurs slightly above or below 
the central line y = 0.5 of the channel. More precisely, let the point of 
contact at the moment of collision have coordinate y = 0.5 + 1. At the 
same time we suppose that the particles have opposite velocity vectors, 
v = (u,v) and — v = (—it, —v), i.e., their total momentum is zero. 

When the particles collide with the opposite walls, their postcolli- 
sional velocity vectors (u',v') and (—u',—v') will again be opposite. 
But the particles collide with the walls at slightly different moments 
of time, and the time interval between their collisions with the walls 
will be 2\£/v\. If, after those collisions with the walls the particles col- 
lide with each other again, the point of contact will have coordinate 
y = 0.5 + £' with \£'\ = \£v'/v\, i.e., 

(9.8) log K| - log M = log(K|/M) = log(sin /(</?)/ sin </?). 

However, if the particles miss each other, then after their second 
collision with the opposite walls their velocities will be again (it, v) and 
(—it, —v). If the particles collide with each other after that, the point 
of contact will be again the distance \£'\ = \£\ from the central line 
y = 0.5. 

Thus the parity issue again arises and the formula (19.81) applies when- 
ever the number of wall collisions between successive interparticle col- 
lisions is odd; otherwise there is no change, \£'\ = \£\. Thus, in the 
notation of (19.71) we have 

m 

(9.9) log|4r| - log |4| ~ ^log(sin/ fc< (^)/ sin 

i=i 

Again we treat </Vs as independent uniformly distributed random 
variables and verify that the mean value of log (sin /(</?)/ sin 93) is neg- 
ative. Proof of the following lemma will be given in Appendix. 

Lemma 9.2. For both k = 0, 1 we have 

1 r 

- I log (sin /fc^)/ sin y?) dtp = fi 2 < 0. 
71 Jo 

Thus the sum in ( 19. 9 j) approaches —00 linearly in m (i.e., linearly in 
time), hence the magnitude of the positional perturbation £ decreases 
exponentially in time. 

This all verifies the stability of the regime S*. Indeed, our perturba- 
tions by 5v account for two directions transversal to S*, and those by £ 



STABLE REGIMES FOR HARD DISKS 



27 



account for one more, hence we took care of all the three codimensions 
in the (reduced) 6- dimensional phase space. 

We note that the classical (unperturbed) system of N = 2 hard balls 
also leaves the manifold S* invariant. The dynamics within S* easily 
reduce to a periodic Lorentz gas with a single particle and infinite hori- 
zon. The periodic Lorentz gas is strongly hyperbolic and ergodic, thus 
the phase points X G S* have one positive and one negative Lyapunov 
exponents, and typical phase trajectories within the manifold S* fill 
it densely. Interestingly, there are no expansion or contraction in any 
transversal direction to S*, i.e., the points X G S* have exactly one 
positive and one negative exponents with respect to the unperturbed 
dynamics in the entire phase space. This makes the manifold S* ex- 
ceptional, as typical phase points are proven [23] to have two positive 
and two negative Lyapunov exponents, cf. Section [2j 

We now describe what happens under our time-reversible pertur- 
bations. First, all the (previously zero) Lyapunov exponents in the 
directions transversal to S* seem to become negative, which instantly 
makes S* an attractor. (Here we refer to Lyapunov exponents of typi- 
cal points X G S*.) Second, the (previously equivalent to the Lorentz 
gas) dynamics within S* is also perturbed, and it seems to retain its 
hyperbolic character and admit an SRB measure. For periodic Lorentz 
gases with infinite horizon under small perturbations the hyperbolicity 
is proven and a (unique) SRB measure is constructed in [12]. Perhaps 
the arguments of [12] could work in our case, too. 

As a result, the SRB measure living on S* seems to attract nearly 
entire phase space. This makes it (the only) physically observable 
measure, as it describes the distribution of typical phase trajectories. 
Therefore, some kind of chaotic behavior does exist in the present case, 
but only after the six-dimensional phase space is reduced to the 3D 
surface S*. 

Lastly, we tried to find similar stable regimes for N > 3 balls with 
time-reversible twists at the walls, and did not observed any - the 
dynamics seem to be totally chaotic. 
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Appendix 



First we prove Lemma I9TT1 We note that f' k (<f) > and J* f' k ((f) dip - 
ir. With these assumptions the integral over (3 can be taken analyti 
cally. It results in 

»27T 



log 



l-[f'M} sin 2 (/3-^) 



d/3 = 4tt log 



Now Jensen's inequality works to show that 



7T 



log 



dip < log 



log 



dip 







Lemma 19.11 is proved. □ 
Now we prove Lemma 19.21 Our arguments apply to both functions 

fo and fi, so we suppress the index and denote them by /. 

Recall that the graph of the function ip = f(<p) is symmetric about 

the line ip = ir — tp\ see Figure [T3J Thus it is convenient to use new 

variables 

__f(<p) + <p , __f(ip)-ip 



and 



w 



2 2 
If we drop a perpendicular from the point (ip,f(ip)), that lies on the 
graph, onto the main diagonal line ip = <p, then its footpoint will be 
(s, s), and its length will be wy/2. The function ip = f{tp) becomes, 
in new variables, w = w(s), and due to the above symmetry we have 
w(n — s) = w(s), i.e., w is an even function with respect to the center 
s = 7r/2. We also note that w(0) = w(ti) = 0, \w'\ < 1, and 

and 



<fi 



s — w 



fif) = s + w, 



therefore 



dip = (1 — w') ds and 



dip 



df ds 
ds dip 



w 



w 



In these new variables the integral in Lemma [9.21 becomes 



(9.10) 
Note that 



M2 



1 

TT 



log 



sins 



w 



sin (s 



w 



(l-w')ds 



because of symmetry w(ir 
We consider a function 



sin(s + w) , 

log — ds 

sm{s — w) 

s) = w(s). 



(9.11) 



F(t) 



sin s + tw) , , 

log — -w ds 

sin s — tw) 
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Figure 13. The graphs of — f(tp) and w = w(s). 

for < t < 1. Note that F(0) = 0. We will show that F'(t) > 0. 
Indeed, 

2F\t) = / (cot (a + wt) + cot (a - wt)) d(w 2 ) 
Jo 

= w 2 (cot(a + wt) + cot(a — wt))| 

r ( i + tw' i - tw' \ 2 , 

+ / r + -t-^t, r w as 

J \sm [s + wt) sm [s — wt) J 

Recall that w(0) = w(jr) = and \w'\ < 1. This makes the middle line 
vanish. The integrand is always positive, which proves that F'(t) > 0. 
Then F(t) > for all t > 0, and 

/i 2 = -±F(l)<0. 
Lemma [9.21 is proved. □ 
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